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Abstract 



A practical method is described for computing the generators of the algebra of first 
integrals that works well for a large family of Hopf-Zero singularity systems. A Lie 
algebra description is presented for all volume-preserving classical normal form of Hopf- 
Zero singularity systems. This is a maximal vector space of classical normal forms with 
first integral; this is whence our approach works. Systems with a non-zero condition 
on their quadratic part are considered. The algebra of all first integrals for any such 
system has a unique (modulo scalar multiplication) generator. The infinite level volume- 
preserving parametric normal forms of any non-degenerate perturbation within the Lie 
algebra of any such system is computed, where it can have rich dynamics. The associated 
unique generator of the algebra of first integrals are computed. The symmetry group of the 
infinite level normal forms are also discussed. Some necessary formulas are derived and 
implemented on appropriately modified Flossier and generalized Kuramoto-Sivashinsky 
equations to demonstrate the applicability of our theoretical results. 
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1 Introduction 

A Lie algebraic structure of a family of vector fields is important in the sense that the family is 
invariant under a group of permissible transformations generated by the Lie algebra. Thus, cer- 
tain dynamical properties are preserved within the family. Therefore, in the study of dynamical 
systems, it is fundamentally useful to recognize nonlinear families of vector fields that have Lie 
algebra structure. Once this is accomplished, the transformation group makes a classification 
for the vector fields within the family. This classification can be computed through the infinite 
level normal form theory. This gives rise to the study of singular differential systems which has 
significantly contributed to the bifurcation and stability analysis of many systems with com- 
plex dynamics. A good practical knowledge on the Lie algebraic structure of any singularity is 
necessary for an efficient infinite level normal form computation. There are many results on the 
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infinite level normal forms of planar singularities. However, there are considerably less research 
results on the infinite level normal forms of three-dimensional state space singularities, where 
the system is already on its center manifold, see e.g., [3 El [TQl [29], [30]. The reason is that the 
Lie algebraic structure generated by three dimensional state space is more complicated than 
the Lie algebraic structure generated by planar state space. Findings on Lie subalgebras of 
a singular family of vector fields contribute to an understanding of the dynamical properties 
invariant under the transformation groups, generated by the Lie subalgebras. In other words, 
it is a fruitful tool for finding dominant subfamilies of singular vector fields. This also has 
an application for the study of systems with symmetry; e.g., Zg-symmetry, reversible, volume- 
preserving, conservative and dissipative systems. We introduce two versal Lie subalgebras of 
the Lie algebra generated by Hopf-Zero singularity. The family of vector fields from these 
subalgebras is associated with volume-preserving conservative systems. 

The original motivation for this paper is the transformation of conservative Hopf-Zero sin- 
gular systems into the infinite level conservative normal form systems. Conservative systems 
have many applications in real life problems and Hamiltonian systems are among the most 
prominent examples of conservative systems. There are substantial contributions on Hamilto- 
nian systems in the literature. However, the Hamiltonian structure requires an even dimen- 
sionality of the state space. Thus, an eminent alternative for conservative Hopf-Zero singular 
systems on central manifold, i.e., a three dimensional state space, is to consider vector fields 
with a first integral, also see [22J for relevant details. A normal form computation may not 
destroy certain symmetric structures (e.g., volume-preserving, Hamiltonian) when the trans- 
formation group is generated by an appropriate Lie algebra of symmetric (volume-preserving, 
Hamiltonian) vector fields. However, the set of all classical normal Hopf-Zero singularities with 
a first integral is not closed under the Lie bracket and also is not a vector space. Therefore, 
the idea of this paper is not to consider the set of all vector fields with a first integral. Instead, 
we consider a maximal vector space (that is also maximal Lie algebra) of such vector fields. A 
noble property of this maximal Lie algebra is that it comprises of all volume-preserving Hopf- 
Zero classical normal forms. A second notable property comes from our upcoming results, that 
is, any Hopf-Zero singularity can be transformed into a sum of two normal form vector fields, 
one with a first integral from this algebra, and the other without any first integral, see [T31 E3] • 
This provides a normal form decomposition of vector fields into conservative and nonconserva- 
tive vector fields which has many outstanding applications in the global analysis of nonlinear 
dynamical systems such as Melnikov theory, e.g., see [14, 26-28, 31]. A good understanding of 
both conservative and nonconservative families is enlightening for the dynamics study of the 
decomposed normal form. 

A practical method for finding the generators of the algebra generated by first integrals 
of a given differential system is an indispensable subject. Our approach displays such fruitful 
method that works well for a large family of Hopf-Zero singular systems. The idea is to 
transform the system into the classical normal form and check if it is a volume-preserving 
normal form system. Then, our formulas derive the generators for the algebra of its first 
integrals. Therefore, this method is applicable for all Hopf-Zero singular systems that their 
classical normal form is divergent free. 
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Normal form computation is an essentially renowned tool for revealing hidden symmetries 
of the differential system and then, any further normalization must preserve the disclosed sym- 
metries. Normal form computation method given in [13] may not work well for the family of 
volume-preserving conservative systems introduced in this paper. This is because the hypernor- 
malization approach in [13] does not necessarily preserve the volume-preserving conservative 
structures. Thus, in our approach a vector field is first transformed into the classical normal 
form and then is decomposed into a sum of a volume-preserving conservative vector field and a 
nonconservative vector field from another Lie algebra that is constructed in [14] . Based on this 
observation, an appropriate hypernormalization approach is applied. This paper deals with 
the cases when the nonconservative part is zero. When the conservative part is zero, one must 
follow the method introduced in [T3J . Finally for the cases of both nonzero conservative and 
nonconservative parts, the normal form computation in [13] is implemented. 

In this paper we are concerned with computing the simplest (infinite level) normal form of 
the volume-preserving systems given by 

-_J2(k + l-l)a{x l+1 (y 2 + z 2 ) k - 1 , 

~- ~ E ^a l k x l y(y 2 + z 2 f~ l + £ b™x m z(y 2 + z 2 ) n ~ m , (1.1) 

-- - £ ^a l k x l z(y 2 + z 2 f- 1 - £ b™x m y(y 2 + z 2 ) n ~ m , 



where -1 < I < k, < k, < m < n, (x,y, z) G M 3 , a{,b™ G R, a^ 1 ^ 0,a° = 0, and b° = 1. 
This system is the classical normal form of a large family of Hopf-Zero singularity systems with 
a first integral. The formal function 

f{x,y,z) := ^a l k x l+1 (y 2 + z 2 f~ l+1 (where - 1 < I < k and I + k > 1) (1.2) 

is a first integral for the system ( 11. lj) . The idea is to analyze the dynamics from the infinite 
level normal form and its first integral rather than working with the more complicated system 
and first integral ( ll.11 - ll.2p . 

The rest of this paper is organized as follows. Section[2]discusses the notations and algebraic 
structures associated with the volume-preserving classical Hopf-Zero normal form with a first 
integral. We obtain the infinite level volume-preserving normal forms and associated first 
integrals in Section [3] Furthermore, a technique is demonstrated in which it enables one 
to change certain coefficients into non-algebraic numbers in order to facilitate elimination of 
certain terms. Then after some normal form computation, one changes them back into the 
desired numbers, see Remark 13.21 In Section HJ under some technical conditions, we prove that 
there exist invertible transformations sending any non-degenerate perturbation of the system 
( II. ip into a truncated infinite level parametric normal form (in cylindrical coordinates) 

x = p 2 ± x p+1 + EI =1 x^ +1 m + J2k= P +i UkX k+ \ 

p =T^^ P -E[ =1 ^ i ^P^-pEf =P+ i M ^^, (i-3) 

o -i + EC+i^-^i + E^/W, 

for some r, s,p, q > 1, N — 1 > rij > — 1, and N — 1 > m 8 > 1 (see Equations (l4.5H4.6j) ). where 
«fc = for k = p — 1 mod 2(p + 1), and j3i = for I = — 1 mod 2(p + 1) and for I = p + q 
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mod 2(p + 1). Furthermore, the formal function 

1 r N 

/(*, P ) = p 2 (^ 2 ± X ^ + £ *-+ v> + E «^ fc+1 ) 

1=1 fc=p 

is the unique generator (modulo scalar multiplication) for the algebra of first integrals for the 
infinite level parametric normal form system ( 11.31) . see Theorem 14.31 In Section |5] we provide 
the necessary formulas for a Hopf-Zero singular system whose cubic-truncated first level normal 
form is governed by Equation (II. ip . Using these formulas, we first apply our approach to a 
modified Rossler equation and then to a symmetric system chosen from a commonly used family 
of generalized Kuramoto-Sivashinsky equation. 

2 Volume-preserving conservative vector fields 

In this section, we provide the necessary algebraic structures for the infinite level volume- 
preserving normal form computation of Hopf-Zero singularity that have first integral, see [151 
24"] for terminology and original ideas. We denote 



n : = ^^r( (t+ i-o.|-^|-^|), (-!<«<*> (») 



J^„.2 , _2\fc-J_ ® JL(„2 i Jl\ k ~ l „ " 



el := sv + ^) fc -'*^ - *V + * 2 ) y^ (o<i<k). 

The vector fields F l k and Q l k can also be represented in cylindrical coordinates, but many first 
integrals for Fj. only appear in terms of x, y, z, see Lemma \2. 11 For aO/ a^ 1 E R, let 

J? :=sveay{e° + a^F - l + J2*irt + J2 b n e n\Kn>l,alb™ eR}, (2.2) 

where — 1 < I < k and < m < n. These notations provide a tool to use a similar approach to 
the method developed by Baider and Sanders [3-4]. The space Jzf is a Lie algebra by the Lie 
bracket [v,w] = vw — wv for any v, w G Jf. Let 

1 := span{6[! + ^ b l k <d[ \ < I < k, 1 < k, b[ G R}. 

This is equivalent to the space of phase components in cylindrical coordinates. The set of all 
formal first integrals for v G 5£ is a subalgebra of formal power series in terms of x, y, z. Since 
the algebra of formal power series is a Noetherian ring, it is finitely generated. We denote by 
(pi,P2, ■ ■ ■ ,Pk) the algebra generated by pi,p%, ■ ■ ■ ,Pk G R[[x, y, z}}, unless otherwise is stated. 

Lemma 2.1. For any I and k, let ^ l k be the algebra of formal first integrals for F l k . Then we 
have: 

• The algebra of first integrals for any 7^ v G % is 

(x,y 2 + z 2 ). 
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• The algebra &\ is generated by monomials that their x-degree is (I + 1) and their (y, z)- 
degree is 2{k + 1 — 1), i.e., 

^ = (x l+1 y 2l z 2{k+1 - l) - 21 | i = 0, 1, . . . , k + 1 - /). 

• For any ^ a 6 M, the algebra 

(x l+1 (y 2 + z 2 ) k+1 ' 1 ) 
is the algebra of first integrals for 0{j + aF k . 

Proof. Let g be a formal first integral for 0™. Then, z-^-g — y-j^g = and g = g(x, y 2 + z 2 ) G 
(x, y 2 + z 2 ). The second part is a straightforward computation. Since every formal first integral 
for Gq + aF l k must be a first integral for 0[| {e.g., see [22j Proposition 3]), any first integral for 
6g + aF l k belongs to (x, y 2 + z 2 ) fl &\. Thus, the proof is complete. □ 

The above lemma suggests the following proposition. 

Proposition 2.2. Let v = 0g + o^Fq 1 + £ a l k F l k + £ 6^6^ G JSf, w/iere a^ 1 ^ 0. T/ien, t/iere 
exists a unique formal first integral f (modulo scalar multiplications) such that the algebra of 
first integrals for v is (f) . 

Proof. Define 

/:= / (v(tx,t*y, th),(y 2 + z 2 ,-2xy,-2xz))dt, 



where (-, •) denotes the usual inner product on M. 3 . Then, / = a J Q 1 (y 2 + z 2 ) 2 + Yl a l k x l+1 (y 2 + 
z 2 ) k+1 ~ l . The formal function / is a formal first integral for v because v(f) = 0. By [22| 
Proposition 3] any first integral for v lies in (x, y 2 + z 2 ) and is also a first integral for w := 
a^F^ 1 + ^2 a l k F k . Let g(x, y 2 + z 2 ) be a first integral for w and g = gN + - • • , where g N denotes 
the nonzero homogenous polynomial component of g with the least degree. Since w(g) = 0, 
the least degree of the expansion for w(g) must be zero, i.e., F ~ 1 (p , 7v) = 0. By Lemma 12.1 1 
gN G ((y 2 + z 2 ) 2 ) and there exist natural number k and real number a\ G R such that A^ = 4fci 
and g^ = a\{y 2 + z 2 ) 2kl . Now let g± :— g — f kl . Then, g\ is also a first integral for w whose 
the least degree nonzero homogenous monomial has a degree of strictly greater than N. An 
induction argument proves that there exist real numbers Oj G R and natural numbers k{ (for 
i = 2, 3, . . .) such that g% '■= g — Y^j=i a jf ki - The proof is complete since gi is convergent to 
zero with respect to the filtration topology. □ 

Remark 2.3. Any v G 5£ has a first integral and is volume-preserving, i.e., div(f) = 0. 
However, ££ is not the set of all vector fields with a first integral. Indeed, any function of the 
form g(y 2 + z 2 ) is a first integral for 

. (I + 1) ,/ d d d \ ., . ,, i d 
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for any I G N. However, the vector fields wi — Fj and 

K*a = _2i»v + *')(2.| + »| + *|) 

do not have any first integral and are not volume-preserving, see [HI Thm 2.3]. This implies 
that the set of all vector fields with a first integral is not closed under the Lie bracket and it 
is not a vector space. Indeed, any classical normal form with Hopf-Zero singularity is uniquely 
decomposed into a volume-preserving vector field (with a first integral) from 5£ and a non- 
volume-preserving vector field (without any first integral), see [13]. This implies that ££ is the 
set of all volume-preserving Hopf-Zero classical normal forms. 

The following lemma portrays the structure constants involved in this paper. 

Lemma 2.4. The following holds true. 

[F l k ,F™] = (( m + i)(k + 2)-(l + l)(n + 2))F%Z, for -I <l <k,-l <m<n, 
[FlQ™] = (m(H2)- n (!+l))e£, for - 1 < I < k,0 < m < n, 

[9^,9™] = 0, for 0<l<k,0<m<n. 

The space X is a nontrivial Lie ideal (and a trivial Lie subalgebra) for Jzf . The above lemma 
implies that the quotient Lie algebra |r is Lie-isomorphic to a proper Lie subalgebra for the 
one-degree of freedom Hamiltonian vector fields, see [U Thm 3.7] and [13] for more detailed 
discussion. 

Notation 2.5. Throughout this paper, we use Pochhammer fc-symbol notation, that is, 

(a) k :=a(a + b) (a + 26) . . . (a + (k - 1)6) , 

for any natural number k and real number b. 

3 The infinite level normal forms 

In this section, we obtain the simplest normal form of Hopf-Zero systems given by Equation 
(II. ip . Let v = YlT=o Vk b e a Lie-graded expansion of v G £6 '. Define d k ' 1 : Jzffc — > j£f& by 
d k ' l (Y k ) = [Y k , Vq] and then inductively define the maps d k,n : S£ k x kerrf fc, " _1 — )■ Jzf/t by 

n-l 

d k ' n (Y k , F fc _!, . . . , F fc _ n+ i) := J2^ Yk '^ u *l' for k - n - 

i=0 

For k < n, define d k,n := d n,n . A normal form style is a rule on how to choose a unique 
complement space C k,n so that Imd k,n © C k,n = ££ k . Then, for a given vector field v, a graded 
Lie algebra structure and a normal form style, there exists a formal invertible transformation 
that transforms v into its n-th (infinite or simplest) level normal form w where w = YlT=o w k 
and w k € C k,n (w k G C k,k ) for all k, see [15] for details. 
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Lemma 3.1. There exists an invertible transformation transforming v^ into its second level 
normal form 

oo oo 

v^ = ° + aoF- + £ a^F k k + £ ^ej, (3.1) 

fc=i fc=i 

where a = a^ 1 . 

Proof. Define a grading function by S(Fl) := 5(6^.) := fc. Then, the result is deduced from 
[F - 1 ,^] = 2(/ + l)F^ 1 and[F " 1 ,ey = 2/ei- 1 . ' D 

Note that by a linear change of state variable, we can remove 0[j from the system, see [HI 
Thm 4.1] and J23J Lemma 5.3.6]. Assume there exist a k ^ (for k > 1) and let 

p:=min{A;|4 2) ^ 0, A; > l}. (3.2) 

Let a := |» a p := a P anc ^ ^p := I-^o" 1 + a pFp- Define a grading structure by 

5(FJ C ) := p(k - I) + k and 5(Q[) := p(k - I) + k + p + 1. 

Remark 3.2. For any a, b £ F, a, 6 > 0, through linear changes of variables 

1 aba sign.(a a p ) _l_ a&a sign(a a P K-!- v v , v 

t := - — ( )v+ 1 t,x := ( )p+ 1 X,y := Y, and 2 := Z, 

ba>Q a p a p 

we can transform SC V into 3£ p := 6F ~ + asign(a p ao)FP. Consequently, the coefficient of F _1 
and F? can be arbitrarily chosen in the normal form computation. Thus without the loss of 
generality, we can choose a := |. Consequently, this enables the coefficient change of a p into 
a non-algebraic number. We shall use this in Theorem 13.51 in order to simplify the system as 
desired, and then again we change it into q p = ±1. Asa theoretical result, this may not violate 
the principles of normal form theory, but this needs more attention when it is implemented in 
a computer program. In fact irrational numbers are treated like rational numbers in computers 
because of computers' round off errors. This, however, does not hamper our results. Indeed 
for implementation of the results on any computer, one needs to truncate the system up to 
a certain degree and thus, one only needs to choose a p to be distanced from the roots of a 
finite number of polynomials. For a case of Bogdanov-Takens singularity when certain ratio of 
coefficients is non-algebraic the simplest normal form is obtained while the case is still open for 
when it is algebraic, e.g., see [Hj. Hence, this technique is very useful wherever it is applicable. 
Indeed, we believe that it can be applied in different other cases, e.g., see 



The following lemma introduces the transformation needed for elimination of the term O™ 
when m < n. 

Lemma 3.3. For natural numbers m and n, there is a 5-homogenous polynomial state solution 
y™ such that 



,n— m— 1 

p.m , ri;m ay 1 _ _J_ m ; '_2 ; H ''V+ 1 of 

^n i Wn ■> <^p\ I \„_ m / . -i\/, , o i „\n—m-l ^pn-pm+n- 



'(vim — n + 1) — n) (m + 1 + (m — n + l)(v+ In 

^ y ' I \ v ' yl ' /p+l fopn—pm+n 

(-2)«- m (m + l)(m + 2+p)£;r " 
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Proof. The proof is straightforward by 

n _^_i (_iy+i ap i ( m + i + ( m _ n 

i=0 



y 



in 
n 



l + 2)(p+l))! 



2p+2 n m+ip+i+l 



2 i+1 (m+ l)(m + p+2) 



^e 



p+1 



n+ip 



u 



(3.3) 



Lemma 3.4. T/ie (p + 1)-£/j /eve/ normal form of v associated with Equation II 1.1}) is 

v {p+1) = e° + \f« 1 + * v fi + £ o% +x) F k k + £ /t +1) ©£, 

where the first summation is over k ^ p — 1 mod 2(p + 1) anc/ t/ie second summation is over 
kj^-i mod2(p+l). 

Proof. Since F k ,Q k G kerad F i, we follow Baider and Sanders [5] and define ^ := ad(Fg) o 
ad(<^,). Then, 

Sf(fJ) = 4(/ + l)(/-fc-2)^ + 2 ap (fc + l-0(p-^ + (fc-0(P + l))^ +1 , 
^(6™) = 4m(m-n-l)Q™ + 2a p {m-n)(m(p + 2)-n(j) + l))e p p X™ +1 . 

Let {i 7 ^ | — 1 < / < k} be an ordered basis for £f , where its ordering is partially defined by 
F l k -< F™ and Q l k -< Q™ if k < n. So, the matrix representation of £f is lower triangular. Thus, 
for any natural number k there exist ^-homogenous polynomial vector fields 



■yk 

r° 

fc,p 



2k a p m (2k+l)™ 2 mp+m _i 
m=0 2 m (m+l)! 2k+mp— V 

Z7 

2fc+mp+2 



E. 

Efc /fc\ a p fc ~ m 7 7 m(p+l)-l 
m=0 VmJ 2 m 2fc+rr 

y^fc Qp m (fc)^ Q m(p+l) 

Z-^m=0 



*(.F, 



*- x ) 



•1 mod2(p+l' 



2 m m! fc+mp 



5(X£)=p mod2(p+l), 
5(^ = mod p+1, 



so that 



ker(^) = span{ J^ 1 , * p fc , 7^° | fc G N}. 



(3.4) 



(3.5) 



On the other hand, 



["Op, /fc iP J 
P>-Mo J 



o, 

fc(p+l)a P *(fc)-2 Q fc(p+i)+p 



2 k k\ 



6 



fc(p+i)+p' 



I - 2k)a p (2k~l)_ 2 ^(p+ij+p-i 



2 2k (2k)\ 



-F, 



2fc(p+l)+p-l 



Thus, F™, 0™ G Imd m ' p+1 for any m = p — 1 mod 2(p + 1) and n 



m 



n + p + 1 . This completes the proof. 



■1 mod 2{p + 1), where 

D 



Theorem 3.5. There exist invertible transformations (including linear time rescaling) sending 
v given by Equation II 1.1}) into the (p + q + 2)-th level normal form system 



x 

y 



(y 2 + z 2 )±x^ 1 + xP +1 ZZi^ +P 



x 



(p+i) 



oo (k+p+l)a k+p ^ k 



zg( x ) =f ^Y+xPy - x p y J2T=i 

-yg( X ) T ^XP Z ~ X?Z J2Zl {k+P+ l )ak+ * x\ 



(3.6) 
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for p, q > 1, and 



DC 



g(x) :=l + x q Y^P q +kX k 



fc=0 



where a k + p = /or k = — 1 mod 2(p + 1), and /3 g+fc = /or k = p mod 2(p + 1) and /or 
A; = — (g + 1) mod 2{p + 1). Furthermore, the (p + q + 2)-th level normal form system /J3.6]) is 
the infinite level normal form. Let 

1 oo 

f(x,y,z) : = ( 2/ 2 + ^)(-( 2/ 2 + ^)±^ +1 + x p+1 ^« fe+p x fc ). 

fc=i 

and & denote for the algebra of first integrals ofv^°°\ Then, J£~ = (/) and the symmetry group 
of v^°°' is generated by JPQq. 

Proof. Assume pf ^ for some i > 1. Then, define 

q := min {i | (3\ p+l) ^0,i> l}. 

For any k > we have A^ +1 = X]m!=o a p m ( m*) ^zk+mp > an< ^ by Lemma 13.31 there exists a 
state solution Y such that 

[x p k+ \ei} + [Y^ p ] 

= ^ ^ + J) (2A: + 2 - m)<— +i (g + 2(m - fc)(p + l))^"* 1 2k{p+1)+p+q 



By Remark 13.21 without the loss of generality, we may assume that a p is not an algebraic 
number and thus, O™ E lm d m+p+1,p+q+2 , where m = p + q mod 2(p + 1). For k = 21, we have 

n, P = (V 2 + **) W + z 2 ) + «^ P+1 )'(^ - V^) e ker ad^. 
Then, 7^ p is extended to a symmetry for ^, + X]fc= P +i a kFk> n EN, i.e., 

Tl n := (^ + ^((^ + ^) + i:« fc ^)'(,|--^)Gkerad^ +EILpQfcF ,. 
Since T$ n is convergent to /'Oq with respect to filtration topology, the proof is complete. □ 

4 Parametric normal forms 

In this section, we deal with the parametric normal form of a multiple-parametric perturba- 
tion of the system (II. ip . Roughly speaking, this section provides an infinite level parametric 
normal form for the miniversal unfolding for the system ( II .ip , also see [HJ [15l [16j [2U [25] and 
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Remark 14.21 Since nonlinear time rescaling destroys the symmetry of the system (i.e., volume- 
preserving), we do not use nonlinear time rescaling. However, parametric time rescaling is 
permitted when it does not depend on the state variables. Consider a parametric vector field 

w(x, y, z, /i) := J2 <nH^ + E fo in e X> ( 41 ) 

where n = (rii,n 2 , . . . ,n m ), rij GNU {0} for j = 1,2,... ,n m , J2T=i n j > 0, fc > —1, k > I > 
— 1, j > i > 0, 6[} 7^ OjO^q 7^ 0, az\ = 0, and a^ = 0. We call any parametric vector 
field w given in Equation (14. II) . a multi-parametric deformation for v when t> = t/;(x, y, z, 0). 
By a similar argument to Lemma 13.11 and a parametric time rescaling, any multi-parametric 
deformation w for v associated with Equation (jl.ip can be transformed to the second level 
extended partial parametric normal form 

- (2) = E C®> n + E «oXv n + E «S^v n + E 0$e^", 

where /3q q = 1, ao, n = §, fc > — 1, and j > 1. Following a similar argument to [TU Thm 4.1] and 
[23} Lemma 5.3.6], we consider a parametric change of state variable [x, y, z] := exp((Y-^- — 

Z^)th(n)) [X, Y, Z], where h(ji) = 1 + £ Po*L» n and [x, y, z] := x£ + y± + z£ denotes the 
new variables and [X, Y, Z] stands for the old variables. Once all linear terms are omitted, 
using a parametric time rescaling, we can transform w^ into 



w^ 



\k' + E «£^v n + E © e X> ( 4 - 2 ) 

(2) ~ (2) 

where k > — 1 and j ' > 1 . Let Qj^O for some k > 1 , and /3 ■ o ^0 for some j ' > 1 and denote 

p := min {/c | a k ' ^ 0, fc > l} and g := min {j | /3-g 7^ 0,j > l}. 

We define the grading function by 

5(F l k /j n ):=p(k-l) + k+(p+q + 3)\n\ and 6(Q l k /j, n ) := p(k-l) + k+p+l + (p+q + 3)\n\. (4.3) 

Lemma 4.1. There exist invertible changes of variables that they transform w^ 2 ' given by 
Equation fl^.jg| ) into the (p + q + 2)-i/i /eve/ extended partial parametric normal form 



w (p+q+2) .__ 



l -F l + J2 «M^V + E ^ne> n , (4-4) 



where k > —l,j > 1, a^o = /or a// k < p, a p0 = =tl> a?fc,n — /or fe = p — 1 mod 2(p + 1) 
and (3k, a = /or both k = — 1 mod 2(jo + 1) and fc = — (g + 1) mod 2(p + 1). 

Proof. The proof is straightforward by similar arguments in the proofs of Lemma 13.41 and 
Theorem 13.51 Note that the number q must be updated in the (p + l)-th level normal form by 
g:=min{j|/^ o +1 V0,J>l}. D 
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By omitting terms of degree higher (standard degree of polynomials) than or equal to n+ 1 
of a vector field v, we obtain its n-degree truncated (n-jet) vector field and denote it by J n {y). 

For any natural number N, let k := L gf+iT -l anc ^ I '■= N — 2k{p + 1) (where \_a\ denotes 
the integer part of the real number a) and r := k(2p + 1) + / + 2. Denote rij (j = 2, . . . , r) for 
all natural numbers in which 

— 1 < n.j < N — 1 and for all tij > p we have rij ^ p — 1 mod 2(p + 1). (4-5) 

Denote m^ (k = 1, . . . s) for all natural numbers 

1 < TUk < N where m^ ^ —1 mod 2{p + 1) and m^ ^ p + q mod 2(p + 1). (4-6) 

Obviously n\ := — 1, ri2 := 0, 77.3 := 1 and m\ := 1. Now for the parametric normal form w <yP+q+2 ' > 
given in Lemma [4. II we denote 

A N := (aij), where a^ = ajg$ for j = 1, . . . , r, and a^ = P%£?} tlH for j = r + 1, . . . , r + s. 

The truncated normal form J Ar + 1 (iy(p+9+ 2 )) i s called a non- degenerate perturbation when 

rank(Ajv) = r + s. (4.7) 

Remark 4.2. When N = p, a transformation of the form a; = X + e (for an appropriate formal 
function e = e(/x)) can be applied on Equation (I4.4H in order to simplify E a p -x^F^Z\ £t n . 
Then, this reduces the number of parameters left at the truncated parametric normal form. 
Thereby, one needs a weaker rank condition than the rank condition (14 .7p in order that a p- 
degree truncated parametric normal form system would be parametric generic. Indeed, the 
sequence rij skips p — 1 and we have r := k{2p + 1) + / + 1. Then, rank(^4 p ) := r + s. These 
types of transformations have been implemented by Murdock and Malonza [25]. However, we 
shall not use them in this paper. 

The condition (14. 7p guarantees that a reparametrization sends the parametric vector field 
jN+i ^ w (p+q+2)^ j n £ ^g ^j\j _|_ ^.degree truncated infinite level parametric normal form 

1 r N N r+s 

jAH-i^M) ;= -F - 1 ±F; + ^F n > 4 + £ an F: + J2fc® k k+ £ e ^^- ( 4 - 8 ) 

j=l n=p+l fc=q i=r+l 

This implies the following theorem. 

Theorem 4.3. There exist formal invertible changes of state variables and parametric time 
rescalings that send any non- degenerate deformation of the system U.l\) into the (p + q + 2)-th 
level parametric normal form w^ p+q+2 \ and then a reparametrization sends J N+l (w <yP+q+2 " > ) into 
the (N + l)-degree truncated infinite level parametric normal form 

= (y 2 + z 2 ) ± X p +1 + El=i x n ' +l ^ + x p+1 Ef=T a k +P x k , 

= zg(x) T ^x p y - EI =1 ^^VlH ~ x p y Ef=7 ^ k+p+ ^ x \ (4.9) 

= -yg( X ) T i£±±l X P Z - JJ i=1 ^^ - xP Z J2k=l P (fc+P+ 2 1)Qfc+P X fc , 
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for 

r+s N+l-q 

g{x) := 1 + ^2 x mi - r ^ + x q ^ &+*&*, P,q>l, 

i=r+l fc=0 

where ak = for k = p — 1 mod 2(p+l), and /3i = for 1 = — 1 mod 2(p+l) and for 1= p + q 
mod 2(p + 1). Let & be the algebra of first integral for w^°°\ Then, 



N-p 

a k+p x k 



1 r N—p 

J? = ^ + Z *)(Ltf + Z l)± X P+l + Y,X ni+1 H i + X* +l Y, 

i=l fe=l 

The parametric symmetry group of w^°°' is generated by J^Oq. 

Proof. The proof is complete by applying a linear change of state variables [x,y, z] := exp((F^ — 
Z-^-)t)[X, Y, Z] on w <yP+q+2 \ where X, Y, Z are the new variables and x, y, z show the old vari- 
ables. This puts the omitted linear part Gq back into the system. □ 

5 Examples 

In this section, we first derive several relations between coefficients that are enough so that 
the cubic-degree truncated classical normal form system would be in .if . We also obtain some 
useful formulas for the infinite level normal forms of a generic cubic classical normal form 
system v G Jzf . Then, we apply these relations to obtain a one-parameter family of modified 
Rossler and a one-parameter family of generalized Kuramoto-Sivashinsky equations. Then, we 
apply our results on these systems. The system of Hopf-Zero singularity can have very rich 
dynamics. Some major contributions have been made in the literature, see e.g., [9, 17, 19-21]. 
However, there are many degenerate or symmetric cases that their dynamics have not yet been 
investigated. The dynamics study of our examples is performed in [12] . 
Consider a differential system governed by 

\ / a,,-fc \ 

x l y 3 z k . (5.1) 



Let i/ 1 ) denote the cubic-truncated classical normal form of the system ( 15. ip . Then, the fol- 
lowing three relations are the necessary conditions for v^ to be in «£?. 

= aoo2 + &oii + cioi, (5-2) 

= 2a 00 3 + 3&oi2 + 3cio2 - 6002 (3&no - 4oioi + 6c 20 o) - c 02(4a ii - 66 02 o - 3cn ), 
= 6020 + 16a 02 i + 26 2 io + 2c 2 oo + 2c 02 o + 8ai 2 - 166 02 o«ioi + 16c 02 oaon 
— 7con(a2oo — ^020) — 16gho(cioi + 7&oii) — ^ioi(7o2oo — 15ao2o) 
-26 2 oo(^no + 2c 200 + 26 2o + 80101) + 2c 20 o(ciio + 80110 + 8a n). 
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Then, 

t,(D = 6° + o^ 1 + a\Fl + <M + a^F 2 2 + b\Q\ + 6°e? + 6^, (5.3) 

where a = \a 02 o + ^a 20 o, 
i 1 l X 

a i — „Oon + -Cioi, 

Ol — 7777^020 + 2&210 + 2C200 + 2Co20 + Coil(a200 — O020) + a.no(feoii — Cioi)) 

111 

+ T77 C 20o(CllO — 2&20o) + 7777^101 ( a 200 — Oo2o) — 7776200(6110 + 26o2o), 

lb 32 lo 

«2 = n( c 102 + 6012 + C(X)2( C 110 — 20oil + 26o2o) + 6o02(2aioi — 2C200 — ^110)) • 

These equations are derived from the fundamentally useful formulas given in PQ, where explicit 
formulas for 6}, 6°, and b\ are also given. 

Proposition 5.1. Consider the vector field 

v := 9° + a F^ + <M + a\Fx + 4*% + & i i + & ? ? + b l®l 
where ao,a\ 7^ 0. Then, the infinite level normal form is governed by 

U Q U 1 /nil t|l ■« oP ^ r,! 1 <l /FT 1 I i|I ■*' 

^ V2a}|a «il 2 ° a i a i 4V2a{|a a{| 2 

where ±1 represents sign(ao a i)- 
Proof. By state change of variables, we can transform the system into 

v := aoFo" 1 + a\Fl + aJF 2 + f*F° + b\Q 1A + 4 ^ + ^ 9 a 8 . 

4a 4a 

The linear changes of variables 

* := (ggg^l)l T> x := (gggg^)ix,y := y, and z := Z 

a a{ 2a{ 

transform this system into the desired form, where a second linear changes of variables is needed 
to add 0° back into the system. □ 

Example 5.2. Consider a modified Rossler equation governed by 



X - 


= -y-z + dy 2 , 


y = 


= x + ay + ez 3 , 


z = 


= bx — cz + xz + rz 3 . 
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The parameter values d = e = r = gives rise to the well-known Rossler equation and for 
two sets of parameter values, it has Hopf-Zero singularity at origin; one is for parameter values 
a = c, b = l,a 2 < 2, while the other has simple dynamics, see [TJ. Hence, we choose 

7-17 2-1 1 r\ S ^ 2 

a := c, o := 1, d := a — 1, and e := —2a -r H — -ra H — -a, 

17 17 17 



where 

-891a 5 + 510a 2 - 1316a 3 + 1058a + 510a 7 - 170a 4 - 340 



r :- 



and - a/2 < a < \f2. 



15(a 2 -2)(15a 4 -2a 2 -96) 
Then, the cubic-truncated classical normal form is given by 

V W = 6° + aoF,- 1 + a\Fl + a?F x + a 2 F 2 2 + 6?0? + 6}0} + 6 2 9 2 , (5.4) 

where 

a i _ a 2 15a 7 - 24a 5 - 5a 4 - 39a 3 + 15a 2 + 17a - 10 

a ° '~ ~V2^?' ^ : ~ 27T^7 ' ° 2 : ~ 2.5(15a 4 -2a 2 -96)(2-a 2 ) 3 / 2 ' 

15a 9 + 577a 7 + 60a 6 - 1788a 5 - 420a 4 + 504a 3 + 840a 2 + 240a - 480 



a 



o ._ 



' 40(15a 4 -2a 2 -96)(2-a 2 ) 3 / 2 

4185a 8 - 7671a 6 - 1020a 5 - 12331a 4 + 3060a 3 + 19938a 2 - 2040a - 15840 



h\ : 



-240(a 2 - 2) 2 (15a 4 - 2a 2 - 96) 



2 17235a 8 - 28066a 6 - 2720a 5 - 72356a 4 + 8160a 3 + 96888a 2 - 5440a - 32640 

2 -320(a 2 - 2) 2 (15a 4 - 2a 2 - 96) ' 

and b\ := 3 " a ~ 4 . Then, the infinite level normal form for any arbitrary a^0 (a 2 < 2) is obtained 
through Proposition 15.11 

Example 5.3. Consider a generalized family of Kuramoto-Sivashinsky equation given by 

d 
h t + —G(h,h x ,h xx ) = 0, (5.5) 

where 

A common approach is to consider traveling wave solutions h(x, t) = u(x + at). Then, this gives 
rise to an ordinary differential equation 

au + Auu + u" + m (4) + 2au'u" + 2bu"u {3) + 2cu'u" 2 + c(w')V 3) + 2duu'u" + du 2 u" = 0. 

Integrating this equation (assuming, u = is a solution for all a), we have 

au + 2u 2 + u + it 3 ' + a(u) 2 + b(u") 2 + c(u x ) 2 u xx + du xx u 2 = 0. (5.6) 
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Let a := (this is to get a Hopf-Zero singularity), and v := (Vl,i>2, V3) = (u,u',u"). Then, 
using a linear change of complex variables, we may send the system to 

= -{x + y + zf{2 + d(y + zj) - (y - zf{a - c{y + zj) + b{y + z) 2 , 

= iy + (x + y + z) 2 (l + l(y + z)) + \{y - z) 2 {a - c(y + z)) - \{y + zf , 

= -iz + (x + y + z) 2 (l + f (y + 2)) + i(j/ - ^) 2 (a - c(y + z)) - \{y + ^) 2 , 

where (0,0,0) is an equilibrium with eigenvalues {0, ±i}. The parameter values a = b = c = 
d = leads to the well-known Kuramoto-Sivashinsky equation, where it is a volume-preserving 
system. Chang [6] remarked that the quadratic-truncated classical normal form of Kuramoto- 
Sivashinsky equation is in JC, i.e., a volume-preserving vector field with a first integral. How- 
ever, classical normal forms of volume-preserving Hopf-Zero singular systems are not generally 
a volume-preserving system. Indeed, a higher degree-truncated normal form of Kuramoto- 
Sivashinsky system is not volume-preserving and also does not have a first integral. Thus, we 
consider parameter values 

7 11 

b ■= 3a, c := —(31 + 859a - 906a 2 ), and d := -— c + 10a - 10a 2 . (5.7) 

975 7 

Then the cubic-truncated classical normal form v^ 1 ' and its infinite level normal form v ^ are 
as follows. 

m ^n , x„ 1 176a 2 + 161a - l n , .. 1162a 2 - 1068a - 62 „ 

V W ■= e2 + (a - 2)F ( T 1 + F,° + (1 - a)F} + F 2 

-r V ; 1300 l -r V ; 1 -r 325 2 

/ 14333 5737a 1013a 2 \^ n 1 rt1 1901a 2 - 4689a + 2399 ^ 2 

+ 1 6? + -a9| H e 2 ., 

V 62400 62400 3900 / l 2 x 1300 2 

ft o , 1 p-i . pi . g 01 4. V2(15384a 3 - 67767a 2 + 94215a - 38384) ^ 2 



^"" ' 2(1 -a) ' -20800(a-l)VI(a-2)(a-l)| 

., .76a 2 + 161a - l)(581a 2 - 534a 



_^2(a - 2)(581a 2 - 534a - 31) 2 (176a 2 + 161a - l)(581a 2 - 534a - 31) : , 
~ 325(a-l)VI(a-2)(a-l)f 2+ 1690000(a - l) 2 3 

where v^ := v+ for a G (1,2), and for a < 1 or a > 2 we have i/°°) := v_ . 
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